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ABSTRACT 

Energy  release  by  aluminum  burning  behind  the  blast  wave  front  produced  by  heavily  aluminized  explosives  was 
investigated.  An  aluminum  evaporation/reaction  model  within  the  multi-phase  flow  was  applied.  The  modeled  HE 
includes  a  significant  percentage  of  aluminum  particles,  whose  long-time  afterburning  and  energy  release  must  be 
considered. 

As  a  first  step  in  the  overall  3-D  comprehensive  methodology  development,  the  evaporation  of  small  aluminum 
particles,  with  5,  50,  and  500pm  diameter,  was  investigated  in  a  1-D  code.  The  aluminum  particles  can  react  with 
oxygen,  water,  and  carbon  dioxide.  The  resulting  pressure  profiles  are  different  from  those  obtained  for  the  no¬ 
reaction  case  or  using  a  one-reaction  which  only  consider  the  reaction  with  oxygen.  Finally  this  aluminum  burning 
model  was  incorporated  into  our  in-house  3-D  code,  FEFLO. 


HEADING 

Aluminum  particles  are  often  mixed  into  explosives  and  solid  propellants  to  increase  their  energy 
release.  The  Aluminum  particles  are  capable  of  burning  behind  the  detonation  front  and  their 
combustion  time  under  such  high  pressure/temperature  conditions  is  significantly  longer  than 
detonation  time.  If  such  heavily  aluminized  HEs/propellants  are  cased,  the  physical  mechanisms 
are  even  more  complex.  The  flow  environment  is  significantly  different  from  bare  charge 
detonation  and  afterburn.  As  long  as  the  case  is  intact,  no  external  air  (oxygen)  is  available,  so 
most  of  the  aluminum  particles  cannot  burn  other  than  in  anaerobic  reaction,  which  is  fairly 
limited  for  explosives  with  high  loading  of  aluminum  particles.  Hence,  establishing  a  valid 
energy  release  model  of  aluminum  burning  behind  the  detonation  front  (afterburning)  is  required 
for  the  accurate  modeling  of  cased  aluminized  explosives.  Developing  an  accurate  numerical  A1 
bumig  model  is  not  trivial  as  the  HE  chemical  reactions  are  very  complicated,  and  the  aluminum 
reaction  controlling  mechanisms  behind  the  detonation  front  have  not  been  yet  fully  established. 

Kim  et  al.  and  Balakrishnan  et  al.  have  reported  the  numerical  simulation  of  the  aluminum 
afterburning  behind  TNT  detonations  [1,  2].  They  applied  the  Khasainov’s  empirical  quasi¬ 
steady  law  [3]  for  aluminum  evaporation  in  a  multi-phase  flow  and  simple  chemical  reaction 
model  using  infinite  chemical  reaction  rates  [4,  5]  to  the  afterburning  behind  the  blast  wave. 
Nobel-Abel  EOS  [6]  was  applied  to  the  main  computation  of  the  blast  wave  propagation  and 
TNT  (C7H5N3O6)  was  assumed  to  decompose  to  4  species  such  as  N2,  H20,  CO,  and  C  as 
follows:  C7H5N206  — » 1.5N,  +2.5H20  +  3.5C0  +  3.5C  [7],  However,  practical 

explosives/propellants  are  often  composites  of  HEs,  oxidizers,  and  binders. 
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The  objective  in  this  study  is  to  develop  a  numerical  model  of  aluminum  particle  burning  behind 
the  blast  wave  (afterburning)  of  heavily  aluminized  explosives  and  to  establish  a  reasonable 
numerical  model  without  resorting  to  expensive  CPU  calculations. 


FLOW  SOLVER  AND  CHEMICAL  REACTION  MODELING 


The  governing  equations  are  as  follows: 
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where  the  three  added  source  term  vectors  are  related  to  the  mass  transfer  (Si),  momentum  (S2), 
and  energy  (S3)  between  the  phases.  The  first  row  of  the  governing  equations  denotes  chemical 
reaction  handling  k-th  reactions.  Second  to  6th  rows  are  the  governing  equations  for  gas-phase, 
and  7th  to  12th  rows  are  the  governing  equations  for  solid-phase.  Subscript  g  denotes  gas-phase, 
subscript  5  denotes  solid-phase,  and  subscript  k  denotes  k-th  species  respectively.  <f>s  denotes  the 
solid-phase  volume  fraction,  and  Np  denotes  particle  number  density.  The  tenns  Ug  and  Us 
denote  velocity  vectors  for  gas  and  solid-phase.  The  terms  Tg,  Ts  also  denote  the  temperatures  of 
the  gas  and  solid  respectively.  The  terms:  5,  d,  and  h  denote  drag  force  factor,  average  particle 
diameter,  and  energy  exchange  factor  between  the  phases  respectively.  They  can  be  written  as 
follows: 


S^Jd)CxPe\ug-u\  (3) 

where  Cx  is  the  drag  force  coefficient  obtained  from  a  formula  proposed  by  Hendrson[8]. 


d  =  (60,  /  xNp)in 

(4) 

h  =  6 </>s(NuAg  Id2) 

(5) 

where  Nu  is  a  Nusselt  number  calculated  from  equations  derived  by  Carlson  and  Hoglund[9]. 
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where  Pr  is  the  Prandtl  number  is  defined  as  follows: 
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where  Cpg  is  the  gas-phase  specific  heats  at  constant  pressure,  jug  is  the  gas-phase  viscosity,  and 
Xg  is  thermal  conductivity. 

The  interphase  mass  exchange  factor  Ac  can  be  written  as  follows: 

Ac  =  (3^s/?s  /r)(l  +  0.276VRe)  forT  ^  T  ignition  (930  K) (8) 

where  Re  is  a  relative  Reynolds  number  based  on  the  particle  diameter  and  velocity  difference 
between  the  gas  and  solid  phase,  ris  a  characteristic  time  of  combustion  described  as  follows: 

z  =  Krdl  (9) 

where  Kr  is  a  burning  rate  constant  and  do  is  an  initial  diameter  of  particles. 


We  considered  4  reactions  and  8  species  [Al,  A1203,  AlO,  02,  H20,  C02,  H2,  C]  in  this  study  as 
follows: 


Al  +  ~^02  =  ^ Al203  -810 kJ / mol  for  T  <  Tdecomposition  (=3500  K)  (10) 

Al  +  ~02  =  AlO  +  66.9  LI /mol  for  T  >  Tdecomposition  (=3500  K)  (1  1) 

Al  +  ^H20  =  ^A1203  + 1 H2  -447.5U /mol  (12) 

A!  +  ^C02  =^Al203  +^C-514.9kJ/mol  (13) 

1-D  COMPUTATIONS 

Al-O?  Detonation 

To  test  the  numerical  model,  an  Al-02  detonation  computation  was  performed.  In  this  test  case, 
detonation  propagation  in  a  dilute  mixture  of  5pm  aluminum  particles  that  were  uniformly 
dispersed  in  an  atmosphere  of  pure  oxygen  was  calculated.  The  solid  phase  concentration  was 
1.5e-3  g/cc.  The  burning  rate  constant  Kr  was  set  to  4.0e6  s/m".  The  mesh  resolution  was  2 
mm/cell.  The  other  initial  conditions  were  shown  in  Table  1. 


Table  1:  Initial  conditions  for  Al-02  detonation  computation 


High  pressure  zone  Low  pressure  zone 


Pressure  100  atm  1  atm 

Gas  phase  temperature  2000  K  300  K 

Gas-phase  velocity  0  m/s  0  m/s 

Solid-phase  temperature  2000  K  300  K 

Solid-phase  velocity  0  m/s  0  m/s 


Figure  1  shows  the  calculated  pressure  profdes  of  AI-O2  detonation  wave,  while  Fig.  2  shows  the 
corresponding  detonation  velocity.  The  detonation  velocity  is  unstable  because  the  combustion 
(reaction)  region  separates  and  catches  up  with  the  shock  front  repeatedly  during  the  propagation. 
The  Chapman-Jouguet  (C-J)  detonation  pressure  obtained  using  Cheetah  is  about  37.1  atm  and 
the  C-J  detonation  velocity  is  about  164600  cm/s.  The  experimental  pressure  and  detonation 
velocity  published  by  Strauss  [10]  are  about  32  atm  and  155000-160000  cm/s  correspondingly. 
The  calculated  pressure  behind  the  peak  (von-Neumann  spike)  is  about  36  atm  as  shown  in  Fig.  1. 
The  calculated  detonation  velocity  is  152000  to  165000  cm/s.  These  calculated  detonation 
pressure  and  velocity  agree  well  with  the  C-J  values  and  the  experimental  data. 


Fig.  1  Pressure  profile  of  Al-CL  detonation 
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Fig.  2  Detonation  velocity  histories 


Highly  Aluminized  HE  Detonation  Modeling 

Here  we  modeled  a  stick  of  highly  aluminized  explosive  with  a  length  of  6cm.  The  JWL  EOS 
[11]  was  used  for  detonation  modeling  from  ignition  to  end  of  detonation.  Cheetah  was  used  to 
determine  the  initial  detonation  product  species.  Assuming  the  Aluminum  particles  were  inert, 
Cheetah  calculated  the  species  fractions  at  frozen  temperature  (1800  K)  as  follows: 

Table  2:  Species  fractions  at  frozen  temperature  (1800  K) 


Name 

Phase 

(mol/kg) 

(mol  gas/mol  HE) 

h2o 

Gas 

7.72E+00 

4.55E-01 

n2 

Gas 

4.38E+00 

2.59E-01 

co2 

Gas 

4.09E+00 

2.41E-01 

hcl 

Gas 

2.55E+00 

1.51E-01 

ch4 

Gas 

2.24E+00 

1.32E-01 

CO 

Gas 

1.72E+00 

1.02E-01 

nh3 

Gas 

5.23E-01 

3.08E-02 

h2 

Gas 

3.55E-01 

2.09E-02 

other  gas  products 

Gas 

3.22E-01 

1.83E-02 

liquid  products 

Liquid 

1.22E+01 

7.22E-01 

solid  products 

solid 

2.61E+00 

1.54E-01 

Total  Gas 
Total  Cond. 


2.39E+01 

1.48E+01 


1.41E+00 

8.76E-01 


From  this  Cheetah  results,  the  8  species  used  in  Eqns.  (10)  -  (13)  were  chosen  as  initial 
detonation  products.  The  JWL  parameters  were  also  obtained  from  Cheetah  calculation. 

Reaction  Model  Effects 

Figure  3  shows  pressure  profiles  at  0.15  msec  for  three  cases  where  different  reactions  are 
considered.  The  A1  particle  diameter  is  set  to  5pm.  The  black  line  (Al  inert)  denotes  the  pressure 
profile  where  aluminum  particles  were  inert.  The  green  line  (one-reaction_A15micron)  denotes 
the  case  where  aluminum  particles  were  only  reacted  with  oxygen  as  per  Eq.  (10).  The  magenta 
line  (Four-reactions_A15micron)  denotes  the  case  where  all  reactions  in  Eqns.  (10)  to  (13)  were 
considered.  The  reaction  in  Eq.  (12)  and  (13)  release  less  energy  than  the  energy  released  in  Eq. 
(10)  as  the  reaction  in  Eq.  (11)  is  endothermic.  Hence  the  pressure  for  the  case  modeling  all  four 
reactions  (Eqns.  (10)  to  (13))  is  less  than  that  for  the  case  considering  only  one  reaction  in  Eq. 
(10). 

Figure  4  shows  the  species  density  profiles  at  0.15  msec  when  all  reactions  were  considered.  The 
A1  particles  around  the  surface  of  the  HE  were  reacted  with  Oxygen  in  the  ambient  region  and 
AlO  was  generated.  The  small  amount  of  generated  AlO  was  located  close  to  the  blast  wave 
front  region.  On  the  other  hand,  most  of  the  A1  particles  were  reacted  with  FEO  and  CO2  in  the 
detonation  products  and  AEO3  was  formed.  Hence,  most  of  the  AEO3  existed  behind  the  blast 
wave  front.  The  amount  of  AlO  is  significantly  smaller  than  other  species.  Hence  AlO  density  is 
plotted  separately  in  Fig.  4. 


Fig.  3  Pressure  profiles  at  0. 15  msec.  1-D  A1  burning  behind  the  blast  wave 
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Fig.  4  Species  densities  at  0.15  msec  1-D  A1  burning  behind  the  blast  wave 


Particle  Diameter  Effect 

Figures  5,  6,  and  7  show  pressure,  impulse,  and  A1  solid  phase  concentration  profdes  at  0.15 
msec  for  explosives  that  contained  different  particle  diameter  sizes.  All  reactions  in  Eqns.  (10)  - 
(13)  are  considered  for  this  test.  The  black,  magenta,  and  green  lines  in  Figs.  5-7  show  the  results 
for  5pm,  50pm  500pm  diameters  respectively.  The  light  blue  line  in  Fig.  7  shows  the  A1  solid 
phase  concentration  initial  value.  The  5pm  A1  particles  burned  quickly  and  produced  the  highest 
pressure  values  at  0.15  msec  as  shown  in  Fig.  5.  On  the  other  hand,  the  500pm  A1  particles  did 
not  burn  completely  and  resulted  in  the  lowest  pressure  values.  Figure  8  shows  pressure  and 
impulse  comparison  among  the  different  particle  diameter  cases  when  the  blast  front  reached  100 


cm  from  ignition  point.  The  difference  of  evaporating  and  burning  velocity  caused  the  different 
pressure  and  impulse  profdes  shown  here.  Tables  3  and  4  show  the  pressure  and  impulse  values 
at  ignition  point  (X  =  0).  The  5pm  particles  burned  quickly  and  resulted  in  the  highest  pressure 
and  impulse  value  at  this  location.  Since  the  50pm  particles  and  500pm  particles  burned  slower 
as  they  propagated  behind  the  blast  wave,  they  released  energy  farther  downstream.  The  pressure 
value  of  50pm  case  at  ignition  point  is  higher  than  that  of  5pm  case  at  this  time  because  the  A1 
particles  are  still  burning,  as  shown  in  Figs.  6  and  7.  Never  the  less,  table  4  shows  that  as  all 
wave  propagated  a  distance  of  100cm,  the  complete  release  of  energy  from  the  5  pm  particles 
yields  higher  pressure  and  impulse  than  the  still  burning  50pm  particles,  which,  in  turn  is  higher 
than  the  slower  burning  500pm  particles. 


Fig. 5  Pressure  profdes  at  0.15  msec  for  3  different  particle  diameters  (5,  50,  and  500pm) 
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Fig. 6  Impulse  profiles  at  0.15  msec  for  3  different  particle  diameters  (5,  50,  and  500pm) 


Fig. 7  Aluminum  solid  phase  concentration  profiles  at  0. 15  msec  for  3  different  particle 

diameters  (5,  50,  and  500pm) 


Fig.  8  Pressure  and  Impulse  profiles  when  the  blast  front  reached  100cm  from  the  ignition 
point  for  3  different  particle  diameters  (5,  50,  and  500pm) 


Table  3.  Pressure  and  Impulse  value  at  ignition  point  at  0.15  msec  for  3  different  particle 

diameters  (5,  50,  and  500pm) 


Time  =  0.15  msec 

Pressure  (dyne) 

Impulse  (dyne.s) 

5  pm  particles 

1.81e9 

3.20e6 

50  pm  particles 

2.1 le9 

2.77e6 

500pm  particles 

1.05e9 

2.27e6 

Table  4.  Pressure  and  Impulse  value  at  ignition  point  when  the  blast  front  reached  100cm  from 
the  ignition  point  for  3  different  particle  diameters  (5,  50,  and  500pm) 


Blast  front  location  =  100cm 

Pressure  (dyne) 

Impulse  (dyne.s) 

5  pm  particles 

1.81e9 

3.20e6 

50  pm  particles 

1.68e9 

2.84e6 

500pm  particles 

7.46e8 

2.31e6 

3-D  COMPUTATIONS 

The  aluminum  burning  model,  which  was  tested  in  1-D  code,  was  incorporated  into  our  in-house 
3D  code,  FEFLO.  Figure  9  shows  the  pressure  Gouraud  shading  of  the  small  tube  test.  The  HE 


Pressure  (dyne) 


modeled  here  is  6  cm  thick  with  average  5pm  A1  particle  diameters.  Figure  10  shows  the 
comparison  of  pressure  profiles  between  results  obtained  with  and  without  the  aluminum  burning. 


Fig.  9  Pressure  Gouraud  shading  at  50  psec 


X  (cm) 


Fig.  10  Pressure  profiles  at  50  psec.  (magenta:  with  A1  particles,  Black:  Inert  A1  particles) 


CONCLUDING  REMARKS 


Energy  release  by  aluminum  burning  behind  the  blast  wave  of  heavily  aluminized  explosives 
was  investigated.  An  aluminum  evaporation/reaction  model  within  the  multi-phase  flow  was 
applied.  The  applied  model  was  validated  via  1-D  Al-O?  detonation  wave  computation.  The 
computed  values  are  agreed  well  with  C-J  values. 

As  a  first  step  in  the  overall  3-D  comprehensive  methodology  development,  the  evaporation  of 
small  aluminum  particles,  with  5,  50,  and  500pm  diameter,  was  investigated  in  a  1-D  code.  The 
aluminum  particles  can  react  with  oxygen,  water,  and  carbon  dioxide.  The  modeled  HE  includes 
a  significant  percentage  of  aluminum  particles,  whose  long-time  afterburning  and  energy  release 
must  be  considered. 

The  resulting  pressure  profiles  are  different  from  those  obtained  for  the  no-reaction  case  or  using 
a  one-reaction  which  only  consider  the  reaction  with  oxygen.  The  A1  particle  size  difference  also 
shows  the  different  pressure  profiles.  While  the  5 pm  particles  were  burned  immediately  behind 
the  blast  wave,  the  50/500pm  particles  were  burned  slower  behind  the  blast  wave. 

Finally  the  aluminum  burning  model  was  incorporated  into  our  in-house  3-D  code,  FEFLO.  In 
the  final  presentation,  more  3-D  computed  results  will  be  presented. 
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